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ABSTRACT 


The  main  object  of  this  thesis  is  to  review  the  work 
relating  to  dual  integral  equations.  Most  of  the  methods  for 
solving  dual  integral  equations,  involving  various  special 
functions  as  kernels,  have  been  outlined  and  various  extensions 
of  earlier  works  indicated,  in  brief. 

We  also  consider  triple  integral  equations  and  discrete 
series  analogues.  In  the  last  chapter  we  introduce  and  solve 
quadruple  integral  equations  by  using  operators  of  fractional 


integration. 
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CHAPTER  I 


INTRODUCTION 


Dual  Integral  Equations  have  their  origin  in  mixed 
boundary  value  problems  arising  in  potential  theory.  Such 
problems  with  axial  symmetry  for  an  infinite  plate  involve 
the  solution  of  Laplace's  equation  in  cylindrical  coordinates 
(p,Z)  in  the  region  0  <  p  <  °°,  0<Z<t  under  given 

conditions  on  the  plane  boundaries  Z  ■  0,  Z  »  t.  If  the 
boundary  conditions  are  such  that  either  the  potential  and 
its  normal  derivative,  or  a  linear  combination  of  the  potential 
and  its  derivative  are  specified  over  the  entire  plane  face  of 
the  plate,  an  easy  way  of  solving  it  is  by  the  use  of  the 
Hankel  transform.  The  use  of  this  transform  eliminates 
the  radial  variable  and  the  problem  reduces  to  the  solution  of 
a  simple  differential  equation  in  one  independent  variable  Z 
with  known  conditions  on  Z  =  0,  Z  =  t. 

On  the  other  hand,  if  the  condition  on  the  boundary 
is  a  mixed  one,  i.e.,  if  the  unknown  function  satisfies  one 
integral  equation  over  the  interval  (0,a),  and  a  different 
equation  over  the  interval  (a,°°),  the  solution  by  a  Hankel 
transform  reduces  to  that  of  a  pair  of  dual  integral  equations. 
For  example,  let  p(r,Z)  be  the  potential  of  a  flat  circular 


I 


electrified  disk  of  conducting  material,  with  centre  at  the 
origin,  and  its  axis  along  the  axis  of  Z.  The  potential 
satisfies  the  partial  differential  equation 


2. 


(1.1) 


1  12. 

r  3r 


0 


Let  p  (s,Z)  denote  the  Hankel  transform  of  p(r,Z)  so  that 


CO 

& 

(1.2)  p  (s,Z)  =  /  rp  (r  ,Z)  Jo(rs)dr.,  Z  >  0 

o 


Differentiating  (1.2)  partially  twice  w.r.t.,  Z,  we  get 


(1.3) 


2  *  oo  2 

2-V  -  /  r  M  J  (rs)dr  . 
3Z2  o  3Z2  0 


We  substitute  the  value  of 
and  get 


from  equation 


(1.1) 


in  (1.3) 


(1.4) 


az  2 


-  -  / 


OO  2 

(r  +  ff)J  (rs)dr. 
o  3r2  3r  0 


Using  the  relation 


CO  2  °° 

/  r  JQ(rs)dr  =  -  /  {  Jq (rs)+rsJ^ (rs)  }dr 


3. 


(1.4)  becomes 


2  *  °° 

3  P  _  f 

2  ■> 

3Z  o 


3p 

3r 


rsJ ' (rs)dr 
o 


00 

=  -s  /  p{ J ? (rs)+rsJ"(rs) }dr 
*  o  o 

o 


(1.5)  =  s  f  prJ  (rs)dr. 

o 

o 


We  substitute  (1.2)  in  (1.5)  and  get 


(1.6) 


2  * 

3  p  2  * 

- T  =  s  p 


Solving  (1.6),  we  have 


*  .  ,  N  -sZ  .  _  .  N  sZ 

p  =  A(s) e  t  B(s)e 


since  p  0,  as  Z  -*  °°  we  have  clearly  B(s)  =  0 


(1.7) 


A(s) 


-sZ 

e 


Applying  Hankel's  theorem  to  (1.7),  we  get 


p(r,Z)  =  /  sA(s)e  sZJ  (rs)ds. 
o 


(1.8) 


4. 


The  boundary  conditions  for  a  disk  of  unit  radius  are 


p  =  constant  (Z  =  0,  0  <  r  <  1) ; 

If  =  o  (Z  =  0 ,  r  ■>  1)  . 


Substituting  sA(s)  =  g(s)  ,  it  follows  from  (1.8)  that  g(s) 
should  satisfy 


CO 

(1.9a)  /  g(s)J  (rs)ds  =  h(r)  0  <  r  <  1 

o 

o 


(1.9b) 


00 

/Qsg(s) JQ(rs)ds  =  0 


r  >  1 


In  general,  the  dual  integral  equations  can  be  written 
in  the  form 


(1.10a)  /  w(y)g(y)K(y ,x)dy  =  f(x)  0  <  x  <  a 

o 


00 

(1.10b)  J  g (y)K(y ,x) dy  =  h(x)  x  >  a 

o 


where  K(y,x)  is  called  the  kernel  function,  w(y)  the  weight 
function;  f(x),  h(x)  are  known  function  and  g  is  to  be 


/ 

5. 


determined , 

In  this  thesis  we  shall  review  the  general  history  and 
development  of  the  theory  of  dual  integral  equations  in  addition 
to  giving  some  extension  of  earlier  works. 

In  particular  we  shall  present  in  Chapter  II  a  historical 
survey  of  the  theory  and  methods  as  well  as  indicate  various 
extensions  of  dual  integral  equations.  In  Chapter  III  we  find 
out  a  solution  of  the  set  of  n-dual  integral  equations  with 
Bessel  function  kernels 

00 

(1.11a)  /  D(a)q(a) J  (at)da  «  f(t)  0  <  t  <  1 

o 

oo 

(1 . lib)  /  E(a)c|)  (a)  J  (at)da  =>  g(t)  t  >  1 

o 

where  D(a)  =  ||D  ||  ,  E(a)  =  ||  E  ||  are  non-singular  matrices 

1C  S  3T  S 

of  known  functions;  r,s  =  l,2J0..,n;  f(t),  g(t)  and  <)>(a) 

are  columns  with  n  coordinates.  The  method  used  is  similar 
to  Szefer’s  method,  in  which  case  Szefer  [77]  has  solved  this 
problem  with  g(t)  =  0.  A  solution  for  the  set  of  n-dual  integral 
equations  with  trigonometric  kernels  can  also  be  obtained  by 
applying  the  same  technique.  In  Chapter  IV,  we  define  quadruple 
integral  equations  by 


(1.12a) 


0  <  p  <  a 


J  5  2%(C)Jv(p5)dS  -  F^p) 

0 


(1.12b)  f  s  2BipU)J  (pz)d£  -  G  (p)  a  <  p  <  b 

v  z 

o 


(1.12c)  |  5  2%(£)Jv(p£)d£  “  F3 (p )  b  <  p  <  c 

o 


(1.12d)  /  C  2%(OJv(pOdC  =  G4(p)  p  >  c 

o 


and  solve  them  by  using  fractional  integrals.  The  solution  is 
obtained  in  terms  of  Fredholm's  integral  equation  of  the  second 
kind . 


The  analysis  used,  throughout  this  thesis,  is  formal 
and  no  attempt  is  made  to  justify  the  integration,  differentiation 
or  the  change  of  the  order  of  integration. 


Finally  we  give  a  fairly  extensive  bibliography  on 


CHAPTER  II 


REVIEW 


Most  of  dual  integral  equations  with  which  we  meet  in 
the  solution  of  mixed  boundary  value  problems  are  of  the  type 


(2.1a)  /  w(x)A(x)K(x,y)dx  =  f(y)  y  e  1^ 

o 

00 

(2.1b)  /  A(x)K(x,y)dx  =  g(y)  y  e  1^ 

o 

where  w(x)  is  a  weight  function;  f(y),  g(y)  are  known  functions 
defined  on  1^  and  1^  respectively  where  1^  =  (0,a)  , 

I2  =  (a,°°);  K(x,y)  is  a  kernel  function  defined  over  the  whole 
xy-plane  and  A(x)  is  a  function  to  be  determined. 

It  seems  that  Weber  [91]  was  the  first  to  consider  dual 
integral  equations.  He  formulated  equations  of  this  type  and  derived 
(by  inspection)  the  solution  for  the  case  in  which  w(x)  =  x  \ 

K(x,y)  =  JQ(xy) ,  f(y)  =  1  and  g(y)  =  0.  Berltrami  [4]  gave  a 
direct  solution  for  the  same  values  of  w(x) ,  K(x,y)  and  g(y) 
but  by  taking  f(y)  arbitrary.  No  work  in  this  field  appears  to 
have  been  done  after  that  for  over  half  a  century,  until 
Titchmarsh  [81]  initiated  the  first  systematic  treatment  for  the 


, 

9  . 


8. 

solution  of  dual  integral  equations.  Applying  Parseval's 
formula  for  Mellin  transforms, 


(2.2) 

°°  .  K+i°° 

/  f (x)g(x)dx  -  /  F(s)G(l-s)ds 

o  K-i00 

where  F(s) 

and  G(s)  are  the  Mellin  transforms  of  f(x)  and 

g(x)  respectively,  to  the  left  hand  sides  of  the  equations 


(2.3a) 

00 

/  f(u)J  (pu)du  =  g(p)  0  <  p  <  1 

0 

(2.3b) 

j  f(u)uJ  (pu)du  “0  p  >  1 

0 

the  value  of  F(s)  is  obtained,  from  which  f (x)  can  be 
determined  by  Mellin* s  inversion  formula.  Titchmarsh  [81]  then 
goes  on  to  solve  the  general  pair  of  dual  integral  equations 


(2.4a) 

00 

/  xaA(x)J  (xy)dx  -  f(y)  0<y  <  1 

0 

(2.4b) 

CO 

/  A(x) (xy)dx  -  g(y)  y  >  1 

0 

with  a  >  0  and  g(y)  =  0.  Busb ridge  [8]  considered  conditions 


for  the  validity  of  Titchmarsh  solution  and  proved  that  if 


' 


9. 


0<a<2  or  -2  <  a  <  0  and  -v-l<a--|-<v  +  l,  then 
the  equations  (2.4)  have  one  and  only  one  solution,  A(x)  ,  in 
the  class  of  functions  whose  Mellin  transforms  are  regular  for 
-v  -  1  <  a  <  a  and  0(|t|  ),  for  every  e  >  0,  in  any 

interior  strip. 

Copson  [19]  suggested  that  if  one  starts  with  a  more 
suitable  form  of  potential  function,  a  single  integral  equation 
is  obtained  involving  a  repeated  integral,  and  this  can  be 
solved  by  two  applications  of  the  known  solution  of  Abel’s 
integral  equation. 

Tranter  [83]  considered  the  dual  integral  equations 
(2.4)  with  a  =  +  1,  g(y)  ^  0.  By  applying  Hankel  inversion 
formula  to  (2.4b)  and  substituting  the  value  of  A(x)  so  obtained 
in  equation  (2.4a),  the  solution  was  obtained  in  the  form  of 
a  Schlomilch’s  integral  equation 


(2.5) 


f  2v+l,  2  2  -^  ,  .  .1  vb  1  . 

j  s  (y  -s  )  X(s)ds  *  (y  tt)  y  P(y) 


where 

(2.6) 


A(x)  =  H(x)  +  (^  irx)'*f  sv+3^x(s) Jv+i,(xs)ds 


and 


' 


10. 


(2.7) 


H(x)  =  g(l)J^+1(x)  +  x  /  yg(y) J^(xy)dy 


Gordon  [36]  solved  the  same  pair  of  dual  integral  equations 
for  general  values  of  a  >  -1.  Assuming  A(x)  =  A-^(x)  +  A^(x) 
and  choosing  A^(x)  such  that 


(2.8) 


!  A2(x)J^(xy)dx 
o 


0 


0  <  y  <  1 


g(y)  y  >  i 


to  get 


(2.9a)  j  x°iA1  (x) J  (xy)dx  =  f(y)  -  /  xaA9(x)J  (xy)dx  0  <  y  <  1 

1  V  £  \) 


(2.9b)  f  A, (x)J  (xy)dx  = 

l  v 


0 


y  >  1. 


A2(x)  is  determined  from  (2.8)  and  A^(x)  from  the  pair  (2.9), 
hence  A(x)  is  found. 


An  alternative  approach  is  to  modify  the  two  integrands 
in  such  a  way  that  they  do  in  fact  become  the  same.  For  this 
purpose,  the  following  Sonine’s  integral  [90] 


(2.10) 


J  ,  (Z)  =  — 

y+v+1 


7y+i  l  4.1  9 

- —  /  j  (Pz)pv  V-pVdp 


2Hr(M+i)  o 


V 


11. 


and  the  discontinuous  integral  [90] 


(2.11) 


f  J, (at) J  (bt)t1+M  Adt  = 

J  A  li 

O 


0 


b|l(a2-bV~|,~1 

^  2Xn,"Vr(A-u) 


0  <  a  <  b 


0  <  b  <  a 


where  A  >  y  >  -1,  are  used.  Noble  [50]  considered  a  more  general 
case  with  g(y)  not  necessarily  equal  to  zero  and  -2  <  a  <  0. 

We  take  f(y)  =  0  and  introduce  a  function  G(x)  such  that 


(2.12) 


00 

/  xG(x)J  (yx)dx 
o 


y  >  1 
0  <  y  <  1 . 


By  Hankel’s  inversion  formula, we  have 


(2.13) 


00 

G(x)  =  /  Cg(C) Jy(Cx)d^ . 


The  integral  equations  (2.4)  with  f(y)  =  0  become 


(2.14a)  /  xai A(x)-xG(x) }J^(yx)dx 

o 

00 

=  -  /  A1+aG(A)J  (yA)dA  0  <  y  <  1 
o 


oc 

f  { A(x) -xG(x) } J  (yx)dx  -  0 
o 


(2.14b) 


y  >  1 . 


12. 


The  dual 

integral  equations  (2.14)  can  be  solved  by  known 

methods , 

The  complete  solution  is  obtained  by  adding  this 

solution 

to  the  one  determined  for  the  case  g(y)  -  0. 

A  direct  method  which  is  an  extension  of  that  given  by 
Copson  [19]  was  considered  by  Noble  [51].  The  use  is  made  of  the 
known  solution  of  Abel’s  integral  equation.  Let  us  set  the 
integral  equation  (2.4b)  equal  to  an  unknown  function  x(y) 
for  (0  <  y  <  1) .  Inverting  by  Hankel's  theorem,  we  get 


(2.15) 

1  00 

A(x)  =  x  }  Ax(A)J  (xA)dX  +  x  f  Ag(A)J  (xA)dA* 

o  1 

Substituting  (2.15)  in  (2.4a),  interchanging  the  order  of 
integration,  we  obtain 


(2.16) 

,,  x  2AV  1  d  f1  ,  2  ,  2^a„,  s  , 

*a)  °  -  r(l+%a)r (jja)  dA  I  s(s  -A  >  G(s)ds 

A 

if  -2  <  a  <  0  and  where 


(2.17) 

x  T(-ha)  s  2v  a  1  d  rS  v+1,  2  2N^ar/  NJ 

<s)  =  r . l4j  ,  ,1+  a  ds  /  y  (s  -y  >  £(y)dy 

F ( 1-H^a )  2  o 

00 

-  /  A1_Vg(A) (A2-s2)"1_^adA  (0  <  S  < 

"  1 

(0  <  s  <  1) 


' 


13. 


The  case  for  0  <  a  <  2  is  treated  similarly. 

Sneddon  [59]  gave  an  elementary  solution  of  the  integral 
equations  (2.4)  in  the  case  a  =  +  1,  v  ■  0,  g(y)  =  0.  For 
a  =  -1,  set 

1 

(2.18)  A(x)  =  x  /  p(t)  cos  xtdt. 

o 

Substituting  this  value  of  A(x)  in  (2.4b),  the  integral  equation 
is  automatically  satisfied,  if  we  use  the  integral  [90] 

]  0  t  <  p 

(2.19)  /  J  (Sp)  sin  (£t)d£  =  ( 

o  I  {  2  2.-h 

[_(t  -p  )  p  <  t 

Now  substitute  (2.18)  in  (2.4a),  then  making  use  of  another  well 
known  integral  [90] 


(2.20)  /  J  (p£)  cos  (£t)d£  = 

y  o 
o 

we  obtain  the  value  of  4>(t)  in  terms  of  Schlomilch's  integral 
equation. 

Copson  [20]  gave  an  elegant  solution  in  the  general  case 


0  p  <  t 

,  2  2,-h 

(p  -t  )  t  <  p 


with  v  $  0,  g(y)  -  0  and  0  <  a  <  2,  -2  <  a  <■  0,  by  putting 


' 


14. 


1-  j  1 

(2.21)  A(x)  =  x  '  (xt)dt. 

a 

o  v+  j 

In  the  course  of  treatment,  solution  of  Abel’s  integral  equation 
and  the  integral  (2.11)  have  been  used. 

Lowengrub  and  Sneddon  [46]  extended  Copson’s  method  [20] 
directly.  They  solved  the  problem  with  f(y)  e  0,  g(y)  4  0  and 
added  this  to  the  Copson's  solution  of  the  pair  (2.4)  with 
g(y)  E  0 .  To  this  end,  we  set 

i  a 
1—  -  00 

(2.22)  A(x)  =  x  ^  f  q> ( t ) J  (xt)dt. 

*  1  v+a 

Using  (2.11)  we  see  that  the  equation  (2.4a)  is  satisfied 
by  this  value  of  A(x)  and  substituting  (2.22)  in  (a. 4b),  we 
get  the  value  of  <J>(t),  so  that  A(x)  is  determined. 

We  observe  that  in  the  methods  given  by  Titchmarsh  [81] , 
Tranter  [83],  Copson  [20],  and  Noble  [51],  the  use  of  integral 
transforms  in  solving  dual  integral  equations  is  quite  complicated 
and  laborious.  The  analysis  is  involved  and  long.  Erdelyi  and 
Sneddon  [30]  introduced  the  operational  approach  for  solution. 

The  fractional  integration  operators  I  ,K  derived  by 

n,a  n,a 

Erdelyi  and  Kober  [29]  and  the  modified  operator  of  Hankel 


transforms  were  used.  They  are  defined  by 


/ 


15. 


(2.23) 

In,af(x)  =  T(a)  (x  y)  y  f(y)dy  a>0,  n  > 

(2.24) 

n  00  _ 

Kn,af(x)  =  1(a)  '  (y  x)  y  f^)dy  a  >  0,  n  > 

X 

(2.25) 

00 

sn,af(x)  ■  x_l5a  fg  y'^J2ll+a(2Axy))f(y)dy 

so  that 


(2.26) 

I  S  =  s 

n+a,3  n>a  n,a+g 

(2.27) 

K  S  0  —  S  ,  0  o 

n,a  n+a,3  n,a+B 

Taking  a  s 

=  23  in  the  integral  equations  (2.4),  we  obtain 

(2.28) 

V+B,-Sf  “  I^v+S,-SS^v-3,2SA  “  S*sv-B,3A 

(2.29) 

K%v-B,38  "  *V-3,3SJsv,oA  *  Si$v-3,3A 

If  we  now  define  a  function  h  by 


(2.30) 

VfB.-U*  °n  (0>1) 

h  “  \ 

V-B,B8  on  (1’°° 

then  h  can  be  calculated  and  hence  A(x)  is  determined. 


« — I  |  C\J  «■ — I  |  CM 


16. 


For  the  case,  when  a  »  -1,  v  =  0,  f(y)  =  0,  in  the 
integral  equations  (2.4),  Lowengrub  and  Sneddon  [45]  found  an 
elementary  method  by  setting 


(2.31)  A(x)  =  x  |  <})(t)  cos  (xt)dt 

1 


where  <j>(t)  is  such  that  the  integral  converges,  this  implies 


(2.32) 


lim  cf>(t)  -  0. 

t  -*■  °° 


Using  the  solution  of  Abel's  integral  equation  and  applying  the 
discontinuous  integrals  (2.19),  (2.20),  the  value  of  4> ( t )  is 
obtained  as 


(2.33) 


2  d 

»(t)  = 


yg(y)dy 

(y2-t2)*4 


Jovin  [39]  gave  conditions  under  which  the  integrands 
on  the  left  hand  side  of  equations  (2.4)  become  the  same  and 
are  reduced  to 


(2.34) 


|  xa/2A(x)JuW2(yx)dx  = 


f  spy) 

\g2(y) 


0  <  y  <  1 

y  >  1 


Multiplying  factor  method  was  used  by  Noble  [53]  to 
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solve  the  pair  of  equations  (2.4).  We  multiply  both  sides  of 
the  integral  equation  (2.4a)  by 


(2.35) 


t  ^  V  1  v+1,  2  2.£ 

-7 -  y  (t  -y  r 

2^r(C+l) 


and  use  Sonine’s  integral  (2.10).  For  equation  (2.4b)  we 
multiply  both  sides  by 


(2.36) 


tv-n-l 

2nr(n+i) 


-v+l 

y 


2  2 

(y  -t  ) 


n 


and  use  another  formula  of  Sonine  [90] 


(2.37)  t"n  1J  _(tx) 

v-n-i 


v-n-1 


2nr  (n+D  o  0 


J  [t(s2+t2)^](s2+t2)  ^vs2n+1ds. 


Integrating  the  resulting  equations  and  inverting  the  orders  of 
integration,  we  solve  the  two  equations  by  taking  the  values  of 
n  in  such  a  way  that  we  have  the  same  integrands  on  the  left 
hand  side.  In  the  special  case  for  f(y)  =  0,  Burlak  [6]  solved 
the  dual  integral  equations  (2.4)  by  transforming  the  two 
equations  so  as  to  involve  the  same  integrands.  Hence  by  Hankel's 
inversion  theorem,  A(x)  can  be  determined.  Recently, 

Srivastava  [72]  has  solved  the  integral  equations  (2.4)  for  the 


' 
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case  when  a  is  a  positive  integer  greater  than  unity.  The  use 
is  made  of  Sonine’s  first  integral  (2.10)  and  Weber-Schaf heitlin 
integral 


(2.38) 


00 

!  Jv+2n+l(x)VXr)dX 

o 


l.i.)  L-  p(o’v)  (2r2-l) 
n+1  n 


0 


0  <  r  <  1 

r  >  1 


where  v  >  -1  and  n  is  a  positive  integer.  The  solution 
is  obtained  in  the  form  of  an  integral  equation  involving  Jacobi 
polynomial  as  kernel. 

William  [94]  reduced  the  dual  integral  equations  (2.4) 
to  the  form 


(2.39a)  F (x)  =  £  iK£)J  (x£)d £  =  f(x)  0  <  x  <  1 

J  n 

o 

00 

(2.39b)  G(x)  =  /  ^ (£) Jn(x£) d£  -  x  ag(x)  x  >  1 

o 


where  0  <  a  <  2  or  -2  <  a  <  0. 

k 

We  denote  the  Mellin  transforms  of  F(x) ,  G(x)  by  F  (s) , 

k 

G  (s)  respectively,  then  from  the  convolution  theorem  of 
Mellin  transforms,  we  have 


F  (s)  =  J  (s)iJj  (1-s+a) 
n 


(2.40) 
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(2.41) 

G  (s)  =  J  (s-a)ip  (1-s-i-a)  . 
n 

Making  repeated  use  of  Mellin  transforms  and  the  Hankel’s 
inversion  theorem,  ^(O  is  determined. 

Uptil  now,  we  have  mostly  concentrated  on  the  solution 
of  dual  integral  equations  (2.4).  After  this,  we  will  consider 
dual  integral  equations  in  which  the  weight  function  is  of  more 
general  nature.  Tranter  [84]  considered  the  following  pair  of 
dual  integral  equations 


(2.42a) 

00 

/  G(u)f(u)J^(pu)du  -  g(p)  0  <  p  <  1 

0 

(2.42b) 

00 

/  f(u)J  (pu)du  =0  p  >  1 

0 

Note  that  in 

(2.4)  ,  G(u)  =  ua  . 

Let  us  assume  that  f(u)  has  a  Neumann  series  representation 


(2.43) 

00 

,,  \  1-K  c  ,  ,  , 

f (u)  =  u  )  a  J„  ,  (u) . 

L  m  2mi-K+v) 
o 

Substituting 

(2.43)  in  (2.42b),  we  see  from 

(2.44) 

00 

I  ul"Kj2m+K+v(u)Jhpu)du  “  0  (p  ’  x) 

o 

that  the  equation  is  satisfied.  On  substituting  (2.43)  in 
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(2 . 42a) ,  the  coefficients  a  are  determined,  after  a  lone 

m 


calculation,  in  the  form  of  an  infinite  series  of  linear  equations. 

Cooke  [14]  found  an  integral  which  is  an  analogue  of 
Tranter’s  solution  [84]  given  by 


(2.45) 


o 


where  0  <  Re(K)  <  1,  and  h(t)  is  a  function  to  be  determined. 
Substituting  (2.45)  in  the  pair  of  equations  (2.42),  the 
second  equation  is  satisfied,  whereas  the  first  one  gives  the 
value  of  h(t)  in  the  form  of  a  Fredholm’s  integral  equation  of 
the  second  kind. 


Burlak  [5]  considered  the  dual  integral  equations  of  the 


form 


cc 


(2.46a) 


u  (u  -K  )  iKu)Jy  (xu)du  -  f(x)  0  <_  x  <_  1 


-y-v,  2  2.  a 


K 


00 


(2.46b) 


/  ip(u)  J^.(xu)du  =  g(x)  x  >  1. 


o 


Since  the  equations  are  linear,  we  may  write 


A(x)  =  A^(x)  +  A2(x) 


21. 


where  A^(x)  denotes  the  solution  in  the  case  g  =  0,  the 

solution  in  the  case  f  5  0.  The  essence  of  the  method  is  the 
reduction  of  the  dual  integral  equations  to  a  single  integral 
equation.  The  use  is  made  of  Sonine's  second  integral  [90] 


(2.47) 


“  J  {a(t2+Z2)S 

/  Jp(bt)  ~2~  2W  tU  dt 
o  y  (t+Z  )2 


0 


0  <  a  <  b 


(2  ,2.J$  v-y-1  9  , 

—  ((  7b  }  }  J  {Z(a2-b2)2}  0  <  b  <  a 

V  v  Z  v-y-1 

a 


where  Re(v)  >  Re(y)  >  -1,  and  of  the  known  solutions  of  the  integral 
equations 

(2.48)  /  ip(p)  (x2-p2)^6J^{K(x2-p2)^}dp  »  h(x) 

a 

and 

00 

(2.49)  /  iKp)(p2-xV6I  {K(p2-xV)dp  “  m(x)  . 

x  6 

Dwivedi  [23]  used  the  multiplying  factor  method  to 
solve  the  equations 

°°  _ii  -V 

(2.50a)  r  u  M  J  (xu){ 1+H(u) }^(u)du  “  f(x)  0  £  x  <  d 

K  M 


22  o 


(2.50b)  /  (u2-K.2)  aJ^ (xu) iJj (u) du  ■  g(x)  x  >  d. 

o 

We  mutliply  both  sides  of  (2.50a)  by 


(2.51) 


,  r  /  2  2 Jix  2  2  p/2  y+1 

Jp{s(p  -X  )  }(p  -X  )^  X  , 


integrate  both  sides  from  0  to  p  and  apply  Sonine's  integral 

2  2 

[90].  On  substituting  s  =  iK  and  / (u  -K  )  =  $ ,  we  have 


(2.52)  /  5  W  P{/(£2+K2)}  V+1J(j+p+1(pO[l+H{v'(C2+Kh}]*{»'(S2+K2)}d5 


2  2. 


-y-p-1  f  cf  \  y+lT  r  T7  /  /  2  2NWv/(p  -x  ,P,  r,  , 

=  p  J  f(x)x  J  {K/(p  -x  )}(  - -)Kdx  0  <  p  <  d 

p  K  — 


provided  p  >  -1,  p  >  -1. 


Similarly,  we  multiply  both  sides  of  (2.50b)  by 


(2.53) 


,  / ,  2  2.,,  y,  2  2  .  -,q  -v+1 

J^{Kv(x  -p  )  }  { '  (x  -p  )}4x 


integrate  from  p  to  00  and  apply  Sonine’s  integral  (2.47).  On 

.,2  2. 

substituting  /(u  -K  )  =  £,  we  get 


00 

(2.54)  /  C2a"V"q{/(C2+K2)}"V"1Jv_q_1(p?)^{/(C2i-K2)}d$ 


o 


.v  9  1  t  V+11  Of7/'v2_n2U^ — 2 — qdx 


=  p  4  f  g(x)x  J  (K/(x  -p  )}(- 


K 


p  >  d 
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provided  v  >  q  >  -1. 

The  values  of  p  and  q  are  so  chosen  that  the  left  hand  sides 
of  (2.52)  and  (2.54)  become  the  same.  The  solution  is  then 
obtained  by  applying  Hankel's  inversion  formula. 

Some  authors  have  solved  dual  integral  equations  with 
trigonometric  kernels  which,  in  general,  are  given  by 

00 

(2.55a)  J  £%(£)  (x£)d5  “  f(x)  0  <  x  <  a 

o 

00 

(2.55b)  /  \p(£)  (x£)d£  «  g(x)  x  >  a. 

o 

Sneddon  [61]  considered  the  system  of  equations  (2.55)  with 
a  =  -1,  by  transforming  the  kernel  functions  into  Bessel  functions 
of  order  one-half.  The  resulting  dual  integral  equations  bear  a 
formal  resemblence  to  the  equations  solved  by  Busbridge  [8], 
However,  the  restriction  -v-l<a--j<l+v  is  violated  and 
the  use  of  this  solution  is  not  justified.  This  fact  was  pointed 
out  by  Fredricks  [35]  who  considered  the  equations  (2.55)  with 
g(x)  =0.  By  applying  Weber  discontinuous  integrals  (2.19)  and 
(2.20),  the  solution  is  determined  only  in  the  case  where  the 
principle  of  superposition  holds. 

Sneddon  [60]  considered  the  equations  (2.55)  with 
g(x)  E  0  and  derived  the  solution  of  the  given  pair  of  equations 
by  using  a  procedure  similar  to  Fredricks.  Here  we  make  use  of 


t 


the  Weber  discontinuous  integrals  (2.19),  (2.20)  and  the 
orthogonal  property  of  the  Jacobi  polynomials.  The  solution 
is  given  by 


24. 


(2.56)  iKO  -  +  2  l  qp  J  (5) 

q=l  "  4 

provided  that  the  constants  p^  are  so  chosen  that 

(2.57)  f  (x)  =  p  +  y  p  P(°’  2)  (x2)  . 

o  L,  q  q 
q=l  n  n 

Following  the  method  used  by  Sneddon  [57],  Srivastava  [70] 
considered  the  dual  integral  equations 


2K-1 


(2.58a)  j  y  g(y)  cos  (xy)dy  =  f(x)  (0  <_  x  <  1) 

o 


(2.58b) 


j  g(y)  cos  (xy)dy  =0  (x  >  1) . 


o 


In  this  case  f(x)  in  (2.57)  is  represented  by  the  Jacobi 
series 


(2.59) 


f(x) 


00 

-  I 

n=l 


a  P(K’^}(x2) 
n  n 


and  satisfies  the  condition 
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(2.60)  /  (l-x^)K  2f(x)dx  -  0, 

o 

Tranter  [88]  solved  the  dual  integral  equations  (2.55) 
by  making  use  of  the  well  known  Bessel  function  Integral 
representations . 


(2.61) 


2 


cos  (£x)  _  f  sin  (gx)dx 


r  **  .  Hx 

j  i  2  2M  dX 

o  (r  -x  ) 


r 


,  2  2M 

(x  -r  ) 


For  cosine  function,  we  integrate  (2.55b)  with  respect  to  ’x’, 

2  2  -k 

multiply  the  resulting  equation  by  (x  -r  )  2  and  integrate 

2  2  -k 

between  r ,  00  .  The  first  equation  is  multiplied  by  (r  -x  )  2 
and  integrated  from  0  to  r,  then  we  get 


r 


(2.62) 


JL  f 

2 


4  k(5)J0(rC)dC 


f (x) dx 


- 


J  2 


2M 


o  (r  -x  ) 


0 


( 0  ~  r  <  1) 

r  1 


Application  of  Hankel's  inversion  theorem  then  gives  the  value  of 

<K5). 

In  the  case  of  sine  kernel,  we  first  differentiate 
(2.55a)  and  then  proceed  in  the  same  way. 

Srivastava  [71]  considered  the  general  pair  of  dual 
integral  equations  given  by 


' 


I  I 
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(2.63a)  J  t^ng(t)  cos  (rt)dt  =  f(r)  0  <  r  <  1 

o 


(2.63b)  /  g(t)  cos  (rt)dt  =  F(r)  r  >  1 

o 

where  n  is  a  positive  integer.  Applying  Fourier  inversion 
theorem  and  the  results  from  Ta  Li[79],  we  get,  under  certain 
assumptions  about  the  nature  of  f(r)  and  F(r) ,  the  value 
of  g(t)  as 


(2.64)  t2ng(t)  =  H(t)  +  (-1)"  |  /  s  2nd{s2n+1q(s)}J2n(st) 

o 


where 


(2.65) 


H(t)  =  j  [t 


it  r .  2n-l  „  ,  ,2n  f 


sin  tF(l)+t  J  F(r)  cos  (rt)dr] , 


(2.66) 


q  ( s)  =  J 


f1  P<'r)T2n+l(r/s) 


,  2  2M 

s  (r  -s  ) 


dr . 


(2.67)  P(r)  =  [f  H(t)  cos  (rt)dt-f (r) ] . 

TT  O 

Dwivedi  [21]  solved  all  the  four  possible  cases  of 
(2.63)  with  sine  and  cosine  functions  by  making  use  of  the 
integral  representations  for  Bessel  functions,  given  by 


, 
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,1-v 


(2.68) 


Vur)  ‘  FfehT)r( 


/UxV  r  cos  (ux)  _  (  .  1 

TX  (p  — ? — TV-m  dx>  Re(v)  p 

^  r  o  (r2-X2p  u  2 


.1+V 


(2.69) 


2 _  ,r.v  t  sin  (ux)  _  ,  .  1 

r&-v)r&)  (u}  ^  (x2-r2)^+v  dX’  Re  v  *  2 


Szefer  [76]  took  the  more  general  case  of  the  dual 


integral  equations 


(2.70a) 


f  G(a)$(a)  (ax)da  =  f(x)  0  <  x  <  1 

4  cos 

o 


(2.70b) 


J  <K°0  (ax)  da  =  0 

4  cos 

o 


X  >  1  . 


Taking  a  trial  solution  as 


(2.71) 


<K°0  =  I  g(S)J  (Sa)d£ 


o 


we  see  that  (2.70b)  is  satisfied.  Substituting  (2.71)  in 
(2.70a)  and  putting 


G(a)  =  a[l+H(a) ] 


we  obtain,  on  using  Weber-Schaf heitlin  integral,  the  solution 
in  the  form  of  a  Fredholm’s  integral  equation  of  the  second  kind. 
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Srivastav  [68]  considered  the  daul  integral  equations 
involving  inverse  Mellin  transforms.  They  are  given  by 


(2.72a) 

-  c+100 

t  — s 

0  .  '  ip(s)  tan  asp  ds  =  f.(p)  0  <  p  <  1 

ZTT 1  .  1 

C-I°c 

(2.72b) 

2Tri  /  s^(s)p  ds  =  f  2  (p )  P  ?  1 

c-i00 

(2.73a) 

Ci-i00 

_  .  /  S\p(s)p  ds  =  f  (p)  0  <  p  <  1 

Z7T 1  .  1 

c-i00 

(2.73b) 

^  c+i00  s 

_  .  /  <Jj(s)  tan  asp  ds  =  f0(p)  p  >  1 

ZTT1  '  ,  l 

C-l°° 

(2.74a) 

1  rC+i°°,  /  n  -s,  =  f..(p)  0  <  p  <  1 

2ri±  !  'PC s)p  ds  1 

c-i00 

(2.74b) 

.  c+i00 

t  —  s 

„  „  j  sip(s)  tan  asp  ds=  f0(p)  p  >  1 

ZTTl  y  Z 

C-loo 

(2.75a) 

..  c+i°° 

f  —  s 

- — -  /  si»(s)  tan  asp  ds  =  f,(p)  0  <  p  <  1 

Z7T 1  J  i 

C-I02 

(2.75b) 

c+i°° 

2jri  /  V^(s)p  ds  =  f 2(p)  p  >  1  . 

c-i00 

, 
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For  the  pair  (2.72)  with  f£(p)  ~  we  assume  that  for  0  <  p  <  1 


(2.76) 


1  fC+1“  ,  -s.  I,1  gl(t) 

m  -  .  a*(s>P  ds  "  -p  37  '  ,V ,2  2. 

C-l“  P  V (t  -p  ) 


dt 


then  from  the  inversion  theorem  for  Mellin  transforms , we  get 


(2.77)  s»(s)  -  '  gl(t)tS  ht  Re(s)  >  -1 


Rewrite  the  equation  (2.72a)  as 


(2.78) 


^  c+i“  ^ 

■j~r  /  ^(s)  tan  —  Trsp  Sds  =  F^(p)  0  <  p  <  1 

C-ioo 


where 


(2.79) 


..  c+i00  - 

If  1  —  s 

Fn(p)  =  f.,(p)  +  J  \p( s)(tan  -z  us-tan  as)p  ds. 


1  2tt  i 


C-ioo 


Solving  (2.78)  and  (2.79),  we  get  the  value  of  g^(t)  in  the 
form  of  a  Fredholm's  integral  equation  of  the  second  kind  given 
by 


(2.80) 


gl(t) 


2  d  ^  pfl(p)dp 

7T  dt  "  ,  2  2 M 

o  (t  -p  ) 


f  g(u) 


u 


K^(u,t)du 


(0  <  t  <  1) 
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where 


c+i°° 


(2.81) 


K^(u , t)  = 


27t1 


f 


c-i° 


/ u  %  s  . 

(-)  (tan 


1 

2 


Trs-tan  as)cot  —  Ttsds. 


Similarly,  the  pair  (2.72)  with  f^(p)  =  0  is  solved  and  the 
two  solutions  are  added  up.  Other  pairs  can  be  solved  in  the  same 
way. 


Srivastav  and  Parihar  [69]  gave  a  closed  form  solution 
for  the  dual  integral  equations,  involving  inverse  Mellin  transform. 
The  following  formulas  from  Mellin’ s  inversion  theorem  and  integral 
relations  [28]  are  used.  For  c  >  0 


for  0  <  x  <  1 
for  x  >  1, 

and  that,  for  c  <  h/2 


(2.82) 


1 


2iri 


C-iloo 


C-l00 


hu  -h 


-1,1  s  -s, 
h  B(—  ,  — )x  ds  = 


(2.83) 


1 

2rr  i 


c+i00 

f 

I 

c-i°° 


,-1^,1  1  sN  -s, 

h  B  ("2 ,  y  -  ^-)x  dx 


fo  for  0  <  x  <  1 

(x^-1)  2  for  x  >  1. 


Set  h  =  tt/os  .  We  assume  that  f2(p)  =  0  in  equation  (2.72b) 
and  that  for  0  <  p  <  1 


1 

2tt  i 


c-fi°° 

t 

J 

c-i00 


siKs)p 


3  t  , .  x  , ,  h  h. 

-p— ;  gl(t)(t-p) 


dt . 


(2.84) 
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Using  inversion  theorem  for  Mellin  transforms  and  formula  (2.8 2), 
we  get 


(2.85) 


-1  ,1  S  V 


s-h/2 


^  Cs)  -  h  B(y,  J  t  gx(t)  dt, 


leading  to  the  relation 


(2.86) 


V(s)  tan  as  =  h  1B(y,-|  -  ^)  J  tS  h/2gi(t)dt. 


Substitution  of  (2.86)  in  (2.72a)  gives  the  value  of  g^(t). 

Erdelyi  [25]  solved  equations  of  type  (2.72)  by  using 
fractional  integration  operators.  The  problem  was  characterised 
in  this  way.  Suppose  that  there  are  two  functions  f,  and  f, 
only  partially  known,  f^(x)  being  given  for  0  <  x  <  1  and 
f2(x)  for  x  >  1;  whose  Mellin  transforms  satisfy  the  relation 


(2.87) 


=  1  (l+n*a“s/m)r  (£+s/n) 
g  2 ( ® =  r(l+n-s/m)r (5+8+s/n) 


with  a(0 . n  known,  then  catl  be  ^eterrain€^* 

Erdogan  and  Bahar  [32]  considered  the  following 
system  of  simultaneous  dual  integral  equations 


*  n 


1  l  a, .  (x)f  (x)J  (xy)dx 


h,  (y)  0  <  y  *  1 


(2.88a) 


' 


, 
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CO  n 

(2.88b)  /  l  b  (x) f  (x)  J  (xy)dx  -  g  (y)  y  >  1 

o  j=i  1J  J  1 


1  =  1,2,.. 


,n. 


By  substituting 


n  °° 

(2.89)  ip  (x)  =  l  b  (x)f  (x)  -  x  /  yg.  (y)J  (xy)dy, 

1  j=l  J  1  1  Mi 


(2.90)  (c..(x))  =  C  =  AB  1,  A  =  (a..),  B  =  (b..), 

ij  i 3  ij 


oo  n  0° 

(2.91)  Pi(y)  =  hi(y)  -  /  J  c± . (x)J  (xy)xdx  /  tg±(t)J 

o  j=l  J  Mi  1  U 

the  system  (2.88)  is  reduced  to 


(2.92a) 


00 


r 


o 


n 

l  c  (x)i^  (x)J  (xy)dx  =  p  (y) 
j=l  3  3  yi 


0  <  y  < 


(2.92b)  j  \jj0(x)J  (xy)dx  =  0  y  >  1; 

o  1  yi 


i  =  1,2, 


9 


n. 


We  set 


(2.93) 


3  00 

VX)  =  X  J  ^  AjmJp  .+2nri-(3  .  <x>  ’  3  =  1>2>- 

J  m=o  j  j 


(xt) dt , 

j 


1 


*  •  5  n » 
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By  this  choice  of  \p  „  (x)  ,  (2.92b)  is  automactically  satisfied 

and  substitution  of  (2.93)  in  (2.92a)  gives  the  value  of 
coefficients  A.  „ 

Westmann  [92]  considered  the  following  pairs  of  dual 
integral  equations 


(2.94a)  /  [a^1(5)-h[»2(5)j J^+2($r)d5  =  f^r)  1  <  r  <  » 

o 


(2.94b)  /  [b^1(^)+^2(5) ] C2aJy+2(5r)dC  -  f2(r)  0  <  r  <  1 

o 


and 


(2.95a)  |  [c^1(^)+^2(0  ]  Jv(Cr)d?=  f3(r)  1  <  r  <  °° 

o 


t  9  a 

(2.95b)  /  [^(O+^tOU  Jy(£r)d£  -  f^(r)  0  <  r  <  1. 

o 

Following  Copson  [20]  and  Lowengrub  and  Sneddon  [43],  the  separate 
cases  are  considered  by  taking  f  (r)  =  0  except  one  of  them 
and  then  assuming  the  representation  for  (£)  ,^2 (£)  in  terms 
of  unknown  functions,  which  are  determined  by  actual  substitution 
of  the  assumed  values  of  iK(£)  and  ^2(£)  in  the  system  of 
equations  (2.94),  (2.95). 
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Szefer  [77]  considered  the  set  of  n-dual  integral  equations 

given  by 


(2.96a)  r  D(a)<f>(a)J  (at)da  =  f(t)  0  <  t  <  1 

P 

o 


(2.96b)  /  E(a)<fi(a)J  (at) da  =  0  t  >  1 

o  P 

where  D(a)  ,  E(a)  are  n  x  n  matrices,  j  E (a)  j  ?  0;  (a)  ,  f(t) 

are  column  matrices  of  order  n0  We  put  E(a)cf>(a)  -  ij;(a),  so  that 

$(a)  ■  E  1(a)i|;  (a)  7  |  E  (a)  |  +  0 

D(a)<j>(a)  =  D(a)E  ^(a)i|;(a)  =  F(a)ip(a) 


and  equations  (2.96)  become 

(2.97a)  !  F(a)ijj(a)J  (at)dt  -  f(t)  0  <  t  <  1 

P 


(2.97b) 


f  ^(a)3  (at)dt  =  0 
P 


t  >  1. 


Again,  we  set 


Y(a) 


!  CP+1g(OJp+1(aC)dC  . 


(2.98) 
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On  substituting  (2.98)  in  (2c97b),  it  is  automatically  satisfied 
in  view  of  the  integral  ( 2 0 11)  and  substitution  in  (2.97a) 
gives  the  value  of  g(£)  in  terms  of  a  sec  of  Fredholm's  integral 
equations  of  the  second  kind,  which  can  be  reduced  to  a  single 
Fredholm's  equation  of  that  kind. 

Erdogan  [31]  considered  the  following  system  of  simultaneous 
integral  equations: 


—  /  P(£)  cos  y£d£  =  pQ(y) 


(y  e  L') 


(2.99) 


f  j  Q(£)  sin  y£d£  =  qQ(y) 

o 


(y  e  L') 


(2.100) 


00 

lim  ^  /  [ (a11+a19x5)P(5)-(a13+a1^x5)Q(C)e  X^sin  y£d£  s  gn (y) 


,  n  7T  -  -  '  11  12 

x  ->  +0  o 


(y  £  L) 


lim  ^  f  [  (a  x+a  xC)P(C)-(a23+a2,xOP(?)e  cos  x£d£,  s  g2(y) 
x  *>  +0  o 

(y  e  L) 


where  L  is  che  union  of  nonintersecting  finite  segments 

L  (L=EL  )  ,  on  0  <  y  <  00  and  L'  is  the  complement  of  L.  Taking 
K  K 

into  account  the  symmetry  conditions  and  setting 


00  CO 

P(£)  =  /  p(t)  cos  £tdt,  Q($)  =  /  q(t)  sin  £tdt 
o  o 


(2.101) 


' 
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the  solution  is  obtained  in  che  form  of  simultaneous  singular 
integral  equations. 

Keer  [41]  considered  the  equations  (2.94),  (2.95). 
Assuming  trial  solutions  analogous  to  Westmann,  the  problem  was 
solved  by  making  use  of  the  operators  of  fractional  integration. 

Another  extension  was  given  by  Fox  [34] .  He  considered 
dual  integral  equations  of  much  more  general  nature  involving 
H-f unction  kernels,  given  by 

00  a  .  ,a  . 

(2.102a)  ■  H(ux  L  i:n)f(u)du  -  g(x)  0  <  x  <  1 

1 p . ,a s 
o  11 

°°  ^4 

(2.102b)  r  H(ux  |  1:n)f(u)du  »  h(x)  x  >  1 

o  Mi,ai 

where 


a,, a  a  ,a1  a  ,a  a  ,a 

(2.103)  H(x  |  6  :n)  -  H(,  |  g  a  :g  a  : ' ' '  :8  ,a  ' 

i  i  112  2  nn 


.  -A-  f 

2tt  i  - 


n  F  (a  tsa  ) 

n  1^7 IT - Mx  S(ls- 

1=1  r(Bi-sal> 


The  path  c 
a  >  -  a,/a,, 

l  l 


is  a  straight  line  x  -  a  and  a  is  such  that 
n  n 

2o  y  a.  <  £  (3c-a0).  We  apply  Parseval's 

i=l  1  i=l  11 


formula  for  Mellin  transforms  (2.2)  to  the  pair  of  integral 


37 . 


equations  (2.102)  and  transform  the  denominator  of  H-function 
in  the  first  equation  to  than  of  the  second  and  the  numerator 
of  the  second  to  that  of  the  first.  Finally,  the  inverse 
Mellin  transform  gives  the  value  of  f(u).  Saxena  [57] 
generalized  Fox’s  result  and  using  the  same  technique  obtained 
a  solution  of  the  following  dual  integral  equations 


(2.104a) 


n  m  (i-a  ,t  ) 

,  H  P’  (xuL  \  /i  r  s  ,,  Q  .  )f  (u)du  =  (x) 

J  0  m,2p+m  1  (y^jou)  ,  (l-d^a^)  ,  (1-3^, t^)  1 


f 


0  <  x  <  1 


00  P+n,0  (y  » X p  ) 

(2.104b)  J  H  .  (xuL,  r  w  \  ,)f(u)du 

'  n,2p+n  1  (A£,S  )  >  (v^a  )  ,  (1-p  ,a  ) 


C2(x) 


x  >  1 


where  HP,q 
n,m 


n  r(e.+B.s)  n  (i-a,-A.s)x 
i  f  i-i  1  1  i-i  11 

27Ti  J  m  n 

L  n  r(i-B.-B.s)  n  r(a.+A.s) 

11  11 
i-p+1  i-q+1 


the  contour  C  being  a  straight  line  x  *  a  and  a  is  such 
that  (i)  a  <  -  3^  /B^  ,  (ii)  o  <  (l-a^)/A^  ,  (iii)  a  fulfills 
the  conditions  of  convergence  ([27],  section  (1.19)). 

Buschman  [9]  pointed  out  that  if  we  define 


ln’a,A(x) 


A 


r(a) 


,  A  .  .a-1  Ap-AaTT/  , 
(x  -1)  x  H(x-l) 


(2.105) 
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where  H(x)  is  the  Heaviside  step  function,  then  the  corresponding 
convolution  theorem  for  Erdelyi-Kober  operators  of  fractional 
integration  is 


(2.106)  Inff  (x)  =  (In,a,A*f)(x) 

x 

where  In,a,A(x)  e  L»(0,°°),  f  £  L'(0,»),  a  >  0,  p  >  j  -1. 

f)  Q(, 

We  have  similar  theorem  for  K  ^  (x) .  Kesarwani  [42]  used  this 

x 

idea  to  solve  dual  integral  equations  involving  Meijer's  function. 
By  application  of  Mellin  transforms  and  fractional  integration 
operators,  the  parameters  in  the  kernels  of  the  two  equations 
are  made  equal  and  then  the  inverse  Mellin  transform  gives  the 
solution.  As  the  dual  integral  equations  (2.102)  and  (2,104) 
can  be  reduced  to  dual  convolution  transform  after  exponential 
change  of  variables,  Tanno  [80]  considered  the  following  equations 


OO 

(2.107a)  j  G (x-t) $ (t) dt  -  f(x)  x  >  A 

—oo 


00 

(2.107b)  /  H(x-t)<j)  (t) dt  =  g(x)  x  <  X 

—00 

where  G(t)  and  H(t)  are  generated  by  certain  meromorphic 
functions  and  A  is  a  constant.  In  order  to  transform  the 
equations  (2.107)  with  a  common  kernel,  we  use  the  uniqueness 


, 
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theorem  of  bilateral  transform  given  by 


(2.108) 

G^G  (x)  =  H^H  (x)  ,  for  all  x, 

* 

where  G^G  is  the  dual  convolution  transform  of  G.  This 
procedure  is  justifiable  only  under  certain  conditions. 

Some  authors  have  considered  triple  integral  equations  of 

the  form 


(2.109a) 

00 

j  <t>(u)J  (ru)du  =  f(r)  0  <  r  <  a 

J  V 

0 

(2.109b) 

oo 

J  u^<£(u)J  (ru)du  =  g(r)  a  <  r  <  b 
o 

(2.109c) 

OO 

J  <J>(u)J  (ru)du  =  h(r)  r  >  b 

0 

where  <J>  (u)  is  unknown.  Tranter  [86]  solved  this  problem  for 
the  case  h(r)  E  0  and  obtained  the  value  of  (u)  in  terms 
of  a  series.  Setting 


(2.110) 

OO 

*(u)  =  U’P  t  Cnb+2n-l+p(bu) 
n=l  r 

we  see  that  the  equation  (2.109c)  is  automatically  satisfied 
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by  this  value.  Substitution  of  (2.110)  in  the  first  two  equations 
leads  to  dual  series  and  the  determination  of  the  coefficients  C 

n 

from  these  series  completes  the  solution.  The  convergence  of 
the  integrals  occuring  in  calculation  requires  that 

111 
Re(v)  > -1  when  p  =  vp  ,  Re(v)  >  -  when  p  =  -  —  . 

Cooke  [15]  made  use  of  the  fact  that  the  equations  are 
linear  and  substituted 


(2.111) 


u2p<J>  (u)  (ru)du  =  f^(r) 

=  f2(r) 


0  <  r  <  a 
r  >  b . 


In  this  way,  we  get  two  pairs  of  dual  integral  equations 


(2.112) 


r 


r  „2p 


u  (u) J  (ru)du  =  f^(r) 

■  g(r) 


/  cf>(u)J^(ru)du  =  h(r) 


0  <  r  <  a 
a  <  r  <  b 

r  >  b 


and 


J  4>(u)J  (ru)du  =  f(r) 


0  <  r  <  a 


(2.113) 


CO 

/  u2p(J)(u)J  (ru)du  =  g(r) 


o 


=  f2(r) 


a  <  r  <  b 

r  >  b 


' 
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From  (2.112)  and  (2.113),  we  find  f^(r)  and  f ^ (r) 
by  applying  Noble*s  method  [53].  So  <J>(u)  can  be  determined 
by  Hankel's  inversion  formula.  The  equations  (2.109)  can  also  be 
solved  by  putting 


$(u)  =  (j^Cu)  +  (P2(u) 


and  writing 

,  00 

j  u2pq1(u)  (ru)du  =  g^r)  0  <  r  <  b 

(2.114)  / 

00 

I  u^^(p^(u)J^(ru)du  —  a  <  r  <  b, 

o 

With  these  substitutions ,  we  again  get  two  simultaneous  pairs  of 
integral  equations,  which  can  be  solved  by  known  methods. 

Cooke  [16]  solved  the  integral  equations  (2.109)  by 
transforming  them  into  the  form  in  which  f(r)  E  0  and  h(r)  =  0, 
This  is  always  possible,  since  by  taking 


(2.115a) 


(  f (r)  0  <  r  <  a 


(2.115b) 


00 

f  $(u)J  (ru)du  =  l  £(r)  a  <  r  <  b 

J  v 

o 


(2.115c) 


V  h(r) 


b  <  r  <  °° 


' 
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$(u)  can  be  determined  by  applying  Hankel  inversion  theorem 
and  so  we  can  find  a  function  m(r) ,  which  is  such  that 


(2.116) 

CO 

[  u2p$(u)J  (ru)du  =  m(r) . 

V 

o 

Write  ip(u)  = 

<p (u)  -  $(u)  and  subtract  equations  (2.115a),  (2.116) 

and  (2.115c) 

from  equations  (2.109a),  (2.109b)  and  (2.109c) 

respectively,  then  we  have 


(2.117a) 

00 

/  0(u)J  (ru)du  =0  0  <  r  <  a 

J  V 

o 

(2.117b) 

CO 

[  u2pb(u)J  (ru)du  =  g(r)  -  m(r)  a  <  r  <  b 

V 

o 

(2.117c) 

00 

f  \p( u)J  (ru)du  =  0  r  >  b 

o 

Now  to  solve  the  equations  (2.117),  we  assume  that 


(2.118) 

oo 

f  \p( u)J  (ru)du  =  x(r)  a  <  r  <  b 

o 

Substitution  of  (2.118)  in  (2.117b)  gives  on  inverting  the 
order  of  integration  and  simplification  a  Fredholm's  integral 


» 
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equation  of  the  first  kind. 

Still  another  method  for  solving  the  integral  equations 

(2.109)  was  given  by  Cooke  [17],  The  problem  is  solved  by 

making  use  of  Erdelyi-Kcber  operators  I  ,  K  and  the  modified 

n,a  n,a 

operator  S  of  Hankel  transforms. 

n  ,a 

The  solutions  obtained  by  various  authors  Williams  [97] 
and  Cooke  [16]  are  reobtained  by  making  use  of  the  operational 
approach  and  the  results  agree  with  those  already  got. 

Dwivedi  [22]  considered  the  more  general  triple  integral 
equations  given  by 


00 

(2.119a)  J  ip(u)J  (xu)du  =  f(x)  0  <  x  <  a 

p* 

o 


oo 

(2.119b)  1  uaip(u)J  (xu)du  =  g(x)  a  <  x  ^  b 

o 


(2.119c)  J  ip (u) JA  (xu) du  =»  h(x)  b  <  x  <  00 

o 


Following  Cooke  [15]  ,  we  substitute  tp(u)  =  ip^(u)  +  ip2(u)  where 
ip^(u)  and  are  determined  by  the  equations  (2.119a), 


(2.119c)  and 
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00 


(2.120a) 


!  u%.(u)J  (xu)du  -  g.  (x)  (0  <  x  <  b) 

1  v  I 


o 


oo 


(2 . 120b) 


u%2(u) J  (xu)du  =  g 2 (x) 


(a  <  x  °°) 


o 


This  gives  us  two  pairs  of  dual  integral  equations  which  can  be 
solved  by  the  analysis  used  by  Noble  [53] .  Hence  the  problem 
is  solved. 


The  discrete  analogues  of  dual  and  triple  Integral 


equations  are  dual  and  triple  series  equations.  Such  equations 
are  obtained,  if  we  make  use  of  series  solution  of  Laplace’s 
equation.  Dual  relations  involving  Fourier-Bessel  series,  Dini 
series,  trigonometric  series,  series  of  Jacobi  polynomials  and 
that  of  Laguerre  polynomials  have  been  solved  by  various  authors. 
For  example,  a  pair  of  dual  equations  involving  series  of  Jacobi 
polynomials  is  given  by 


OO 


(2.121a) 


u=o 


OO 


A 


(2.121b) 


n 


n-o 


L  r(S+n+l)r  (ct+n+1/2) 


where  a 


^  ,  0  >  -1  and  P^a’^(cos  0)  denotes  the  Jacobi 


' 


. 


45. 


polynomial.  Srivastav  [66J  considered  this  problem  in  two 
stages.  In  the  first  case,  we  put  G(0)  a  0  and  assume  that, 
when  0^0  <  <p 


(2.122) 


00  A  t 

\  _  ~  n _ _ _  p(otj3)/  a  \ 

L  r (3+n+l)r (a+n+l/2)  n  1  ; 

n=o 


<f> 


h^(u) du 


,  Ian-23  a  d  f 

=  -(cos  — 0)  cosec  0  —  J  77 - - - v 

2  d0  ^  /(cos  0  -  cos  u) 

0 


Using  orthogonality  relation  for  the  Jacobi  polynomials,  we  get 
on  simplification 


/tt 


(2.123) 


A  = 
n 


qn(a ,3) 
172 


f 


4  h1(u)  sin  (y  u) 

o 


1  .  ,2a-L(a-^,g-^) 


n 


(cos  u)du 


where  q  (a,3)  =  (a+g+2n+l)niF (a+3+n+l) . 

Substituting  the  value  of  from  (2.123)  in  (2.121a)  and 

solving  for  h^(u) ,  we  obtain 


h1(u) 


2  ,  .  1  , l-2a ,  1  v 1+23 

=  —  (sm  u)  (cos  —  u) 

TT  2  2 


d  (U  F(0) (sin  ^0)^asin  0  d0 

du  J  / (cos  0  -  cos  u) 

o 


In  the  second  case,  we  put  F(0)  =  0  and  determine  the 
value  of  A^  in  the  same  way.  The  solution  of  the  general 
problem  is  now  obtained  by  merely  adding  the  two  solutions. 


* 


' 
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Recently,  Askey  [3]  has  considered  a  type  of  discrete  dual 
equations  which  he  calls  dual  sequence  equations.  A  typical 
example  is 

1 

(2.124a)  J  (l-x)°f  (x)P<'a,f3^  (x)  (l-x)a(l+x)^dx  *  a  n  =  0,1,..., N 

*  -1  n  n 

(2.124b)  (  f (x)P^a’^ (x) (l-x)a(l+x)^dx  =  b  n  =  N+l,... 

'  -1  n  n 

The  author  mentioned  that  there  are  a  number  of  results,  given 
by  Feldheim  [33],  Al-Salam  [2]  and  others,  for  the  orthogonal 
polynomials  which  contain  Sonine’s  formulas  (2.10)  and  (2.47) 
as  limiting  cases.  Some  of  these  formulas  can  be  used  to  play 
the  same  role  for  dual  series  and  dual  sequence  equations  that  we 
have  seen  the  Sonine  integrals  to  play  for  dual  integral  equations. 

The  solution  of  dual  sequence  equations  involving  Laguerre  and 
Jacobi  polynomials  has  been  obtained  under  certain  conditions. 


' 


CHAPTER  III 


A  SET  OF  n-DUAL  INTEGRAL  EQUATIONS 


In  this  chapter  we  find  a  solution  for  a  given  set  of 


n-dual  integral  equations.  Making  use  of  matrix  notation  and 


applying  the  discontinuous  integral 


(3.1) 


r  J  (at)J  „ (ax)da  = 


00 


V  v-1 

o 


V-1  ,  _V 

X  it 


0 


t  <  X 


t  •>  X 


provided  Re(v)  >  0, 

the  solution  is  obtained  in  terms  of  a  Fredholm's  integral  equation 
of  the  second  type.  We  shall  employ  a  method  similar  to  Szefer  [76], 
It  is  primarily  based  on  the  method  used  by  Copson  [20]  and  later 
on  extended  by  Lowengrub  and  Sneddon  [46], 

1.  Transformation  into  canonical  form.  Let  us  consider  the 
set  of  n-dual  integral  equations  with  Bessel  function  kernels. 


00 


(3.2a) 


f  D(a)c}>(a) J^(at)da  =  f(t)  0  <  t  <  1 


o 


CO 


1  E(a) $ (a) J  (at)da  =  g(t) 

v 


o 


(3.2b) 


t  >  1 


A 


> 
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where  D(a)  —  |j  D  ||  ,  E(a)  =  j|  E  ||  are  non-singular  matrices 

rs  3T  s 

of  known  functions;  r,s  =  l,2,...,n;  f(t),  g(t)  are  column 
vectors  whose  elements  are  known  functions  (with  n  coordinates) ; 
<{>(a)  —  a  column  vector  whose  elements  are  unknown  functions  (with 
n  coordinates)  and  J^(at)  is  the  Bessel  function  of  the  first 
kind  and  of  order  \)  .  Szefer  [77]  has  solved  the  problem  with 
g(t)  =  0.  Here  we  shall  consider  the  set  of  equations 


(3.3a)  /  D(a)cf>(a)J  (at)da  =  0  0  <  t  <  1 

o 

00 

(3.3b)  j  E(a)$(a)J  (ctt)da  «  g(t)  t  >  1 

o 

since  once  the  solution  of  this  set  of  equations  is  known,  the 
solution  of  the  set  (3.2)  can  be  obtained  merely  by  adding 
Szefer's  solution  to  it. 

In  order  to  solve  the  set  of  equations  (3.3),  we  shall 
transform  them  into  a  convenient  form  by  substituting 

(3.4)  D(a)cf)  (a)  =  x(°0  • 

The  relations  (3.4)  give  us  a  column  vector  of  unknown 
functions  (with  n  coordinates)  and  since  D(a)  is  non¬ 
singular  so  that  it  has  an  inverse.  Thus  we  can  write  (3.4)  in 


' 
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the  form 

(3.5)  cf>  (a)  =  D *  i(a)x(ot)  . 


Equations  (3.3),  can  be  rewritten  as 


(3.6a)  f  y(a)3  (at)da  =0  0  <  t  <  1 

V 

o 

oo 

(3.6b)  J  L(a)x(ct)J  (at)da  =  g(t)  t  >  1 

o 

where  the  functional  matrix  L(a)  is  given  by 

L(a)  =  E(a)D  ^ (a) . 

The  set  of  equations  (3.6)  will  be  called  the  set  of  dual 
integral  equations  in  the  canonical  form  and  it  is  this  set  whose 

solution  will  be  sought. 

2,  Solution  of  the  equations  (3.6).  We  seek  the  solution 
of  the  set  of  equations  (3.6)  in  the  form 

OO 

(3.7)  X(°0  “  t V  /  &(x)J  ,(ax)dx  , 

1 

where  &(x)  is  an  unknown  vector.  Substituting  the  value  of 


iojosv  nwon^ou  n-; 
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x(a)  from  (3.7)  in  (3.6a),  we  get 


OO  CO 

f  tV  /  £(x)J  t (ax) J  (at)dxda  0  <  t  <  1. 

J  J  ,  v-1  v 

o  1 


Inverting  the  order  of  integration,  it  becomes 


(3.8) 


OO  oo 

/  tV£(x)dx  ^  J  (at)J  ..(ax)da  0  <  t  <  1. 

4  V  v-1 


Since  x  >  1,  we  observe  that  t  <  x  always;  hence  applying  the 
discontinuous  integral  (3.1),  the  integral  (3„8)  vanishes 
identically  and  consequently,  the  equations  (3.6a)  are  satisfied. 


Now  we  substitute  (3.7)  in  (3.6b),  then  we  have 


00  oo 

(3.9)  ,r  L(a)[tV  J  £(x)J  n  (ax)dx] J  (at)da  =  g(t)  t  >  1. 

J  ,  v-1  v 

o  1 


Assuming  for  the  matrix  L(a) 


(3.10)  L(a)  »  I  +  V(a) 


where  I  is  the  unit  matrix  and 


we  obtain  from  (3.9)  that 


00  CO 

(3.12)  g(t)  =  f  I  [  t  ^  1  &(x)J  (at)J  n(ax)dx]da 

,  v  v-l 

o  1 

00  00 

+  J  V(a)[tv  f  i  (x)  J  (at)J  _^(ax)dx]da  . 
o  1 


Inverting  the  order  of  integration  in  (3.12),  we  have 


00  00 

g(t)  =  /  tV£(x)dx  /  J  (at)J  (ax) da 
y  y  v  v-l 

1  o 

00  00 

+  J  [J  V(a)J  (at)J  -  (ax) da ]  tV&  (x) dx . 

1\)  \)  «L 

O 


Again  applying  (3.1),  we  obtain 


t  OO  00 

(3.13)  g(t)  =  r  xV  1Ji(x)dx  +  tV  J  [J  V(a)J  (at)J  (ax)  da  ]  i  (x)  dx . 

y .  y  .  y  v  v-l 

1  1  o 


Putting 
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(3.14) 


/  ^£(x)dx  =  Q(t) 

1 


and  assuming  that 


(3.15)  lim  Q(t)  =  0 

t  ->  00 


we  have 


(3.16) 


tU  1J’(t)  =  Q’(t)  . 


Substituting  (3.14)  and  (3.16)  in  (3.13),  we  obtain 


CO  00 

(3.17)  Q ( t)  +  tV  |  [J  V(a)J  (at)J  _.(ax)da]x1  VQ’(x)dx  =  g(t) . 

IV  V““X 

o 


Putting 


(3.18)  /  V(a)J  (at)J  (ax)da  =  N(x,t), 

v  v-i 
o 


equations  (3.17)  may  be  written  as 


(3.19) 


Q(t)  +  tv 


oo 


r 

j 

i 


N(x,t)x^  VQ 


(x) dx  =  g(t) . 


Integrating  by  parts  and  making  use  of  the  condition  (3.15),  we 


get 
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(3.20)  Q(t)  -  tV  /  [N(x,t)x^  V] 'Q(x)dx  =  g(t). 

1 


Taking 


(3.21) 


[N(x, t)x1  V] '  =  -  K(x, t) , 


we  finally  obtain 

cc 

(3.22)  Q(t)  +  tV  /  K(x,t)Q(x)dx  ■  g(t). 

1 

This  is  a  set  of  Fredholm's  integral  equations  of  the  second 
ind,  which  can  be  reduced  to  a  single  equation  of  that  type. 
Hence  the  problem  is  formally  solved. 

Note  that  for  n  »  1,  the  problem  reduces  to  Tranter's 


case  [ 84] . 


CHAPTER  IV 


QUADRUPLE  INTEGRAL  EQUATIONS 


In  the  analysis  of  mixed  boundary  value  problems  in 
which  different  conditions  hold  over  two  different  parts  of  the 
same  boundary  can  often  be  easily  reduced  to  the  solution  of 
dual  integral  equations.  Such  equations  involving  Bessel  or 
trigonometric  functions  as  kernels  have  now  amassed  a  considerable 
literature.  In  the  current  decade  Tranter  [87],  Cooke  [15], 

Dwivedi  [22]  and  a  few  other  authors  have  considered  the  problems 
with  different  conditions  over  three  different  parts  of  the  same 
boundary.  The  integral  equations  involved  in  the  solution  of  such 
problems  were  termed  as  triple  integral  equations.  Ln  some  boundary 
value  problems ,  however ,  the  boundary  condition  is  such  that 
different  conditions  hold  over  four  different  parts  of  the  boundary 
and ,  in  such  cases ,  the  integral  equations  involved  can  be  written 
in  the  symmetrical  form 


t  —  2.QL 

(4.1a)  L1(a,p)  =  /  £  ip(£) J^CpOdC  =  F^p)  0  <  p  <  a 

o 


L2(B,p)  e  /  $  2%(£)Jv(pOd£  =  G2(p) 

o 


(4.1b) 


a  <  p  <  b 


/ 


> 
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OC 

p  —  2a 

(4.1c)  L3(a,p)  =  /  £  iKC)Jv(pS)d£  =  F3(p)  b  <  p  <  c 

o 

00 

r  -26 

(4. Id)  L^(6,P)  =  J  K  ^(5)dv(p5)dS  =  G^(p)  p  >  c 

o 

where  F^,  G^,  F^,  and  are  known  functions  of  p  and 

ip(0  is  to  be  determined.  We  shall  call  such  equations  as 
' Quadruple  Integral  Equations ' .  In  this  chapter  we  derive  a 
solution  of  the  integral  equations  (4.1),  in  the  case  when 
=  0  s  G^,  by  using  operators  of  fractional  integration  and 
those  of  Hankel  transforms.  We  will  show  that  this  can  be 
extended  to  the  case  in  which  G^  and  G^  are  not  necessarily 
zero . 


1.  Operators .  We  give  here  a  brief  account  of  the 
definitions  and  properties  of  the  operators  used  in  solving  the 
integral  equations  (4.1).  We  define  operators  bI  and 

3.  r)  ^  ot 

by  the  formulae 

c  n  ,a 


b  2x 

(4.2)  I  f (x)  =  ZX 
a  n  ,a 


-2a-2q  b 


r(a) 


/  <xW) 

a 


2  2.a-l  2p+l 


f (u) du 


a  >  0 


I  f  (x) 
a  n  ,a 


-2n-2a-l  ,  b 

x _ d_  f 

f (1+a)  dx  " a 


,2  2.  a  2p+l 

(x  -u  )  u 


f (u)du 


-1  <  a  <  0 


(4.3) 


I 
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(4.4) 


K  f  (x) 
c  n  ,a 


2x 


2n 


r(a) 


t  2  2  a-1  -2a-2p+l 

(u  -x  )  u  f(u)du 


a 


0 


(4.5)  f(x) 

c  n»a 


x 


2n-l 


F ( 1+a )  dx 


,  2  2. a  -2a-2p+l 

(u  -x  )  u  f 


(u)du 


-1  <  a  <  0 


For  a=  0,  these  are  just  identity  operators.  Note  that  with 

X  00 

these  definitions  I  and  K  are  simply  the  Erdelyi-Kober 

op, a  xp,a 

operators  [61] .  In  these  cases  we  will  drop  the  indices  on  the 

left  and  write  them  as  I  and  K  .  We  also  observe  that 

n,a  p,a 

(4.2),  (4.3)  make  sense  if  b  <  x  and  similarly  (4.4),  (4.5) 
are  defined  only  if  c  >  x. 


The  modified  operator 
is  defined  by 


S 

n  ,a 


of  the  Hankel  transforms 


(4.6) 


Sn,a£(x)  ’  2V“  !  cl'aj2n+a(xt)i(t)dt- 


Sneddon  [61]  has  shown  the  following  relations  between  the 
Erdelyi-Kober  and  Hankel  operators 


(4.7) 


I  .  _S  »  S  ,0 

n+a,S  n,a  n,a+B 


K  S 

n,a 


n+a,3 


S 


n  ,a+3 


(4.8) 
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n+a,g  n,a  n,a+g 

p,a  n,a+3  risOi+3 

provided  the  conditions  for  the  existence  of  the  various 
operators  are  satisfied.  The  inverse  operators  are 

(4.11)  bI~1  *  bI 

a  n,a  a  n+a , -a 

(4.12)  dK_1  =  dK 

c  n,ot  c  n+a,-a 


(4.9) 

(4.10) 


(4.13) 


S  1  -  S 

n  ,  a  p-ra  ,  -  a 


We  give  two  lemmas,  which  define  the  product  of  pairs  of 


operators 


£  x  ~  1 

Lemma  4.1.  Let  I  ,  be  operators  as  given  in  (4.2) 

— -  -  e  n  ,  a  d  n  ,  a - c — — — —  — - 

(4.3)  and  (4.11).  Then 


x T -1  f  r  ,  s  2  sm  Tra  -2r\  ,  2  2,  -a  f  (d  -t  )  t  f(t) 

(4.14)  .1  I  f  (x)  ■  - - -  x  (x  -d  )  - r — - — 7T- - 

dr)»aen,a  t  j  2,1 


x  -  t 


provided  x  >  d  >_  f  >  e 


dt 
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Proof :  We  shall  prove  it  in  two  steps.  First,  for 
0  <  a  < 1  and  then  for  -1  <  a  <  0. 

Case  I.  For  0  <  a  <  1,  we  have  from  the  definitions  of 

jl  ^  and  ^1  that 

d  p  ,a  e  p  ,a 


Pf(x)  =  *1  1  fI  f(x)  =  *1  .  fI  f(x) 

d  p  ,a  e  p  ,a  d  p+a,-a  e  p,a 


-2p-l  x 

x  d 


r"  ,  2  2s-a  2n+2a+l  2y  2a  2n  ,l  ,  2  2,a-l 

fa^y  to  j ,  (x  ~y  5  y  f(a)‘  dy  '  (y  ~u  > 

d  e 


2p+l  .  . 

*  u  f(u)du 


If  f  <_  d,  then  this  integral  is  real.  Changing  the  order  of 
integration,  we  get 


(4.15)  Pf(x)  =  rTaffeFcO  (E  !  “2n+lf(u)du  /  y(x2-y2)  “ (y2-u2)a_1dy 


In  the  inner  integral,  we  put 


(4.16) 


2  2  2  ,2  .  2 

y  =  x  cos  6  4-  d  sm  0 


and  use  equation 


(4.17)  2F1(a,b;c;Z) 


ryp  /  \  tt/2  ,  .  .2b-l  2c-2b-l 

21  (c)  f  (sin  9) _ (cos  8) 


Kb)r(c-b) 


(1-Z  sin20)a 


d0 


o 
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of  Erdelyi  [28],  then  the  inner  integral  in  (4.15)  is  equal  to 


(4.18)  (1-cO  1Z1  a2F1(l-a,l-a;2-a;Z) 

where 


Z 


2  j 

x  -d 


2 

2  5 


x  -u 


dZ  =  2x(d2-u2) 

dx  ,2  22 

(x  -u  ) 


Substituting  (4.18)  in  (4.15) 


Pf  (x) 


sin  na  -2r)-l 

iT(l-a)  X 


J  U2n+1f(u)  ^  [Z1  a2F1(l-o,X-o;2-a;Z)] 


2  2. 

2x(d  -u  ) 

,  2  2.2 

(x  -u  ) 


du 


2  sin  tot 

TT 


f  2n+l,,2  2.  ,  . 

r  U _ (d  ~U  )f(u) 

J  ,  2  2.2 

e  (x  -u  ) 


Z  a2F1(l-a,l-a;l-a;Z)du 


Since 


(4.19) 


2F1(l-a,l-a;l-a;Z)  =  (1-Z) 


a-1 


Pf(x)  = 


2  sin  tot  -2n  .  2  ,2. -a  u2p+1^2_^2^  f(u)du 

it  X  U  -a  )  j  2  2 

e  x  -u 


Case  II.  When  -1  <  a  <  0 
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Pf(x)  =  XI  1  fI  f(x)  -  *1  .  fI  f(x) 

d  n,a  e  n,ot  d  n+a,-a  e  n,0i 


2x  2n  fx  ,  2  2,-ot-l  2n+2a+l  y  2n  2“  1  ,  d  ,£  ,  2  2Na  2n+l,,  ,  , 

fT^T  J  .  (X  -y  )  y  r(l-ta')  dy  d7  I  <y  -O  u  f(u)du. 

Cl  6 


Differentiating  the  inner  integral  under  the  integral  sign,  we 
have 


T34T/  \  2x  2n  fX  ,  2  2v-a-l  2y  ^  .  2  2.a-l  2n+lr, 

Pf(x)  =  T^a)  (x  "y  }  dy  I  (y  "u  >  u  f(u)du 


T  (a) 


Changing  the  order  of  integration,  we  obtain 


(4.20)  Pf(x)  = 


2x  2n  rf  2r)+l  .  ,  rX  .  ,  2  2.-a-l,  2  2.a-l, 

J  u  f(u)du  J  2y(x  -y  )  (y  -u  )  dy 


r(a)r(-a) 


2x 


-2n  f 


x  2  2 


T  (a)T (-a) 


*  ^  ^  ^  2  2 
f  2n+l^/  n  i  r  /X  -y  N-a-l,  2  2s-l1/x  -y  N 

j  u  f(u)du  J  (  2  J  2)  (u  -x  )  d (  2  2^ 

d  y  -u  y  -u 


2x 


-2n 


(-a)r (a)T (-a) 


f  ..2n+ljr ,..N  /„2  j2^-a  (d2-u2)a 


j  u  f (u) (x  -d  ) 


2  2 

x  -u 


du 


2  sin  Tra 
7T 


(x2-d2)  a 


f  2p+l,  2  2  o.  .  , 

t  u _ (d  -u  )  f(u)du 

j  2  2 


Hence  the  lemma  is  completely  proved. 


Lemma  4.2.  Let  \  be  operators  as  given  in  (4.4), 

0  n  ,a  x  rjjOi 

(4.5)  and  (4.12).  Then 


. 
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(4.21) 


dK  1  fK  f (x) 
x  n,a  e  q ,a 


2  sin  ua  2n+2a,t2 
- - -  x  (d  -x  ) 

7t 


0  f  ,  2  A2.a  -2a-2r\-tl  ,  . 
2..  -a  f  (t  -d  )  t _ f  (t) 

J  2  2 

e  t  -x 


dt 


provided  x  <  d  <_  e  <  f . 


The  proof  of  this  lemma  is  similar  to  that  of  Lemma  4.1. 


2.  Extensions .  Let  us  consider  quadruple  - integral  equations 


00 

(4.22a)  f  $  (£x)d£  ■  f (x)  0  <  x  <  a 

v 

o 


oo 

(4.22b)  /  iKOJ^(£x)d£  =  g(x)  a  <  x  <  b 

o 


(4.22c)  J  rY*(5)Ju(5x)de  =  hW  b  <  x  <  c 

o 


oo 

(4.22d)  J  ip(£)J  (£x)d£  -  k(x)  c  <  x  <  °° 

v 

o 


We  determine  a  function  <f>(£)  which  is  such  that 


/  <J>(£)Jv(5x)d£  =  s  (x) 

o 


(4.23a) 


0  <  x  <  a 
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(4.23b) 

00 

/  <f>(5)j  (5x)d£  =  g(x)  a  <  x  c  b 

'V 

o 

(4.23c) 

®  m(x)  b  <  x  <  c 

(4.23d) 

=  k(x)  c  <  x  <  00 

where  s(x) 

and  m(x)  are  any  convenient  functions  of  x. 

(f)(5)  can  be  determined  by  usual  Hankel  inversion  theorem  and  so 


we  can  find 

a  function  r(x)  which  is  such  that 

(4.24) 

00 

J  5  Y<t>  (5)  Jv(5x)d£  =  r (x) 

0 

Substitute  x(£)  ~  ^(5)  -  <f>  (C )  and  subtract  equations  (4.24), 
(4.23b),  (4.24)  and  (4.23d)  from  (4.22a),  (4.22b),  (4.22c) 
and  (4.22d)  respectively  and  we  have 


(4.25a) 

oo 

/  5  YX(S) J^(£x)d£  “  f(x)  -  r (x)  0  <  x  <  a 

0 

(4.25b) 

00 

/  x(OJ  (5x)d5  =  0  a  <  X  <  b 

o 

(4.25c) 

00 

/  5  \(5)Jv(£x)d£  =  h(x)  -  r (x)  b  <  x  <  c 

0 

(4.25c) 


b  <  x  <  c 


(4 . 25d) 


00 


c  <  X  <  00 
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/  x(C)Jv(5x)d4  =  0 
o 


These  equations  can  be  solved  by  the  method  given  in  the  next 
section  and  thus  x(£)  is  determined.  Hence  \p(£)  is  found. 

3.  Solution  of  the  equations  (4.1).  We  transform  the  equations 
(4.1)  into  a  form  to  which  the  operational  theory  is  applicable 
by  substituting 


9a 

(4.26)  i KO  *  w(5);  f(p)  -  (2/p)  ZaF(p)  ; 

g(p)  *  (2/p)2  G(p) 


by  means  of  which  we  get 


(4.27a)  L  (o,p)  s  22op"2“  /  C1  % (6) (pS)dS  =  q(p)  0  <  p  <  a 

o 


(4.27b)  L2(B,p)  5  226p  2B  j  C1  2Bt(5)Jv(PC)d5  =  g2(p)  a  <  p  <  b 

o 


(4.27c)  L3(o,p)  ;  22ap  2“  /  61  2a*(5)Jv(pC)d5  =  f3(p)  b  <  p  <  c 

o 


(4.27d)  L4(fj,p)  5  22Bp  26  ;  s1  2B*(0 J^pOdC  »  g4(p) 

o 


p  >  c 


' 
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Let  1^  denote  the  interval  (0,a),  the  interval 
(a,b),  the  interval  (b,c)  and  the  interval  (c,“). 

For  a  function  f  in  L~  we  shall  write  f,  +  f~  +  f„  +  f ,  , 

2  12  3  4 

where 

f.  =  f  on  I,  and  =0  on  I.  i,j“l,2,3,4;i^j 

Using  S-operator  defined  in  (4.6),  the  integral  equations  (4.27) 
reduce  to  the  form 

(4‘28)  V«,2»*(P)  =  f(p) 


(4-29)  ^-6,2(3*  (p)  =  g(p). 

Here  f^,  g^,  f^  and  g^  are  prescribed  but  g^,  f^,  and 

are  to  be  determined.  We  consider  the  case  in  which  G„  =  0  =  G, 

2  4 

and  take  as  trial  solution 


(4.30)  p(p)  =  S^g^Kp)  . 

Substituting  this  value  of  $(p)  in  (4.28),  (4.29)  and  using 
formulas  (4.9),  (4.10),  we  have 


(4.31) 


f  -  T  e 

x%v+8,a-8 
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(4.32) 


g  ”  ^1  O  O  £  » 

b  J$v-3,3-a 


Also,  we  have 


(4.33) 


£  -  Ij  ^  ,  Q  0f  , 

%v+3,a-3 


(4.34) 


£  "  \v-3,6-ag 


We  proceed  to  determine  £  and  the  subscripts  on  all 
the  operators  will  be  dropped  for  brevity  sake.  All  I’s  will 
be  supposed  to  have  subscripts  ^v+3 ,  a-3  understood  and  all 
K’s  to  have  ^v-3,3-ot. 


Evaluate  (4.33)  on  1^,  then 


(4.35) 


-1 

£,  =  I  f.,  ,  i.e.,  £,  is  determined. 

1  o  1  1 


Taking  (4.34)  on  1^  we  have 


(4 . 36) 


OO—l  00—1 

£  =  K  g  -  K  0  -  0. 

4  p  fo4  p 


Evaluating  (4.32)  on  I  we  get 


g„  =  bK£0  +  ^K£„  +  °°K£ . 
Z  p  Z  b  J  c  4 


i.e. 
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0  =  bK£0  +  +  °°K0 

p  2  b  3  c 


which  gives 

(4.37)  £2  -  -  V15m3 


where  p  <  b  <  b  <  c.  Applying  Lemma  4.2  to  (4.37)  we  obtain 


(4.38)  £2(p) 


2  sin  TT(g-a) 

7T 


v-2a  2  2,a-8 

P  Tb  -p  ) 


rC  (t2-b2)  B-at-V+2a+1«-3  (t) 


dt 


Finally  evaluating  (4.31)  on  ,  we  have 


f  3  " 


I£. 


+  bi 


+  ,  II. 


a  2  b  3 


Simplifying  we  get 


(4.39) 


£  =  P1  bf 

j  b  r3 


P  T  — 1 

b 


c^l  b 


Pi'1 


bI£? 
a  2 


Since  f ^  and  are  known  functions.  The  function 

(4.40)  d(P)  -  ^r1f3(P)  -  ^r1  ®U1(p) 


is  known.  Substituting  (4.40)  in  (4.39)  we  have 


. 
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(4.41)  £0(p)  “  d(p)  "  b1  1  aI^2(p) 


Applying  Lemma  4.1  to  the  second  term  on  the  right  hand  side  of 
(4.41)  and  substituting  the  value  of  Jl ^ (p  )  from  (4.38)  we 
obtain 


„  ,  s  i  /  \  i  2  sin  TT(a-e)  -v-2B  ,  2  ,2sg-a 
£3(p)  =  d(p)  +  - - -  P  (p  -b  ) 


C  ,u2  __2Na-g„v+2g+lr  2  sin  Tr(g-a)  _v-2a,u2  _2va-g 


/  (b--y-)“  My  v  '  "*■[ 

a 


TT 


y  (b  -y  ) 


,2  2  g-a  -v+2a+l  , •  . 

c  (t  -b  )  t  J&~(t)  1 

f  - - -  dt]  — - -  dy. 


2  2 
t  -y 


2  2 
p  -y 


Inverting  the  order  of  integration,  we  get 


(4.42)  *  (p)  =  d(p)  -  JC  {p-v-26(p2_b2)S-a 

TT  b 


,2  2N g-a  -v+2a+l  t  ,,2  2.2(a-g)  2v-2a+2g+l 

(t  -b  )  t  J  (b  -y  )  y 


a 


2  22  2  dy}£3(t)dt- 

(t  -y  )(p  -y  ) 


Putting 


4  sin  TT(q-g) 

2 

TT 


X 
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and  expression  within  the  curly  brackets  in  (4.42)  =  -  K(p,t), 
we  obtain 

c 

(4.43)  &3(p)  -  d(p)  +  X  /  K(p,t)£3(t)dt 

b 

which  is  a  Fredholm’s  integral  equation  of  the  second  kind  and 
can  be  solved  by  known  methods.  The  equations  (4.35),  (4.36), 
(4.38),  and  (4.43)  completely  determine  i  and  hence  our 
problem  is  solved. 
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